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Abstract 


This report focuses on the linear analysis of a plane-parallel flow stability in trans¬ 
verse magnetic field (Hartmann flow) within the convective approximation. We obtain 
and solve equations describing the perturbation growth. We found the perturbation 
modes and their non-excitation conditions. We obtain the equation for the instabil¬ 
ity increment and show that it has an instable root. Also we shown that resulting 


instabilities qualitatively agree with the experimental data. 

Introduction 

Hartmann flow is a steady stream between two fixed infinite parallel planes arising due to 
the pressure drop that occurs in a magnetic field directed orthogonally to planes. We choose 
the z axis to be co-directional with the external magnetic field B 0l and the x axis direct along 
the stream (see Fig. [[]). For such a flow there is an exact solution: 



( 1 ) 



where k\ = B 0 /y/Airpvv m , k 2 = — (1 / pv) (dp/dx). Constants v and u m are kinematic and 
magnetic viscosity respectively, p is a density of the fluid. 

From (HD, © we see that with increasing transverse magnetic fields velocity profile be¬ 
comes flatter. This is due to the Lorentz force acting on different areas differently. Where 
the speed of the current is less than the average speed, the electric current flows along the y 
axis and against it when the speed exceeds the average value. In the first case, the Lorentz 
force accelerates the flow and in the second case it slows the flow down. This leads to a 
flattening of the velocity profile. This flatness is characterized by the Hartmann number: 
Ha = B q / y/An pvv m . 

Stability of the Hartmann flow was first considered in [Tj and there was obtained a di¬ 
mensionless equation similar to the hydrodynamic equations of Orr - Sonnnerfeld: 



( 3 ) 
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Figure 1: Profiles of the velocity field for the Poiseuille flow (dash line) and Hartmann flow 
(solid line). 


where u means an unperturbed velocity, ip - stream function perturbation, a and c - dimen¬ 
sionless wave frequency and phase velocity of the perturbation. The boundary conditions 
are ^(±1) = ^'(±1) = 0. The right side of the equation which takes into account the direct 
effect of the magnetic field on the perturbation is negligible. So we obtain the usual equation 
of Orr - Sommerfeld, but with Gartmon velocity profile, which can be replaced by a close 
one: u — 1 — exp(—Ha( 1 + |z|)). Thus the influence of the magnetic field is taken into 
account only by changing the velocity field and we obtain a (quite expected) result that with 
an increase of the magnetic field stability increases too. For sufficiently large (Ha > 20) 
magnetic field there is a linear dependence of critical Reynolds number on the Hartmann 
number: Re c = 5 • 10 4 Ga. 

In [;2] the stability of the ideal conducting fluid between two coaxial rotating cylinders 
is considered. They write down equations describing an evolution of perturbations and by 
analyzing them (but not solving) derive a stability criterion of the current in the layer between 
cylinders (Velikhov criterion) dfl(r)/dr > 0, where H is the angular velocity of the fluid layer. 
Comparing this criterion with the hydrodynamic Rayleigh criterion d(r 2 Q)/dr > 0 we can 
conclude that the magnetic field destabilizes the flow in a rotating cylinder up to a certain 
value of the magnetic field, and strong magnetic field stabilizes the flow. This is due to the 
fact that in a weak magnetic field electrodynamic forces are already affecting the nature of 
small-scale motions, and the effect of freezing has not yet manifested. 

Also worth noting the work of [3] in which the transition to turbulence due to the insta¬ 
bility Hartmann layer and conditions of turbulence suppression (laminarization) have been 
studied experimentally . It has been found that when the parameter R = Re/Ha > 380 
the flow becomes turbulent. Numerical simulations [4] give approximately the same result. 
Because R is inversely proportional to the magnetic field it may be concluded that a weak 
magnetic field destabilizes the current. 

The question then arises: how important is the role of rotation in the formation of mag- 
netorotational instability? Is there a parallel flow instability and under what conditions they 
are excited? Let’s try to answer these questions. 
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Two-dimensional perturbations 


Assume instabilities are convective, i.e. perturbations that arise at any point do not have 
enough time to develop and are carried over beyond the real pipe. But since the magnetic held 
has (because of embeddedness) an inhibitory effect, for a feasibility of such an assumption 
the magnetic held should be small. 

The magnetic held can be compared to the magnetic viscosity by choosing multipliers to 
make the dimensions to match. The same result can be achieved by requiring the Alfven 
speed Bq/ y/Anp to be much smaller than some characteristic velocity of the fluid. We can 
construct this value from dimensional parameters of the liquid in three ways: z/ m /d, u/5 and 
y/vv m l5. Since the embeddedness is affected only by the magnetic viscosity, it is logical to 
choose the hrst way. Thus we obtain: 


BnS 


V / 47rpz/„ 


< 1. 


(4) 


In this approximation we suppose that the perturbation does not evolve and moves along 
the main how, ie d/dx = 0. Then we investigate stability of the system of equations: 
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V x (V x B) + i/ m V 2 B, 


(5) 

( 6 ) 


V-V = 0 , 


(7) 


V • B = 0, (8) 

assuming that the main how obeys equations (jT]), (J2]) and all instabilities are convective. 

We make transformations B— ;• B + b, V— )• V + v, P —> P + ip, where V = (V (z), 0, 0), 
where B = (B(z),0, Bq) is known. Then we leave only linear terms of perturbations. Then, 
since the movement is infinite in time t and coordinate y, we assume that the perturbation 
is periodic in terms of these variables: f(y,z,t ) — > f(z)exp(i r yt — iky). I.e. in each layer 
dz there propagates a plane wave. The boundary conditions for perturbations correspond to 
adhesion and impermeability conditions; the wall is non-magnetic, so perturbations on them 
tend to zero: 


b(±<5) = v(±<5) = (^(±5) = 0. (9) 

Since there is no held sources, it is possible to introduce the vector potential, and as 
d/dx = 0, then y and z rotor components comprise only one term. Denoting then b x = b 
and v x = v we write: 


, da da. 

h = (h 'Tz^ 

dq dq 


( 6 , ^,ika), 

(v,^,ikq). 


( 10 ) 

( 11 ) 
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Due to these transformations equations (JTJ) and (J 8 ]) disappear, and y and z components 
of the equation ([ 6 ]) are identical. Then the system of equations describing the amplitude of 
a plane wave perturbation has the form: 


. B 0 db B 2 d 2 v 

I'yv + ikV q — - - - ik- — a + uk v — u—— = 0, 

47T p dz Anp dz 2 

dv d 2 b 

i'fb — ikV'a — B 0 -—b ikB’q + v m k 2 b — z/ m -—^ = 0, 


dz 


dz 2 


( 12 ) 

(13) 
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— ikN 
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where stand out two quantities: 


A + ikM = 0. 
dz 


M 


i'yq + uk 2 q — v 


d 2 q 
~dz 2’ 


(14) 

(15) 


(16) 
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(p Bb 

-^ 1 —• 

p Anp 


(17) 


For N the boundary conditions are obviously zero: N(±5) = 0. However, another con¬ 
dition is needed to determine the possible wave numbers k. Note that if we make a replace¬ 
ment in the expression for M: z —> —z (or k — y — k ) it does not change. So we can say that 
M{8) = M(—S). This statement can also be proven from similarity of the expression for M 
with the heat equation (in x - space). Then M acts as a heat source. And since the walls 
are identical, they will generate perturbations in the same manner. 

Then from (|T 6 ]h (fI71) and the boundary conditions which are set above we obtain the 
eigenvalues of the wave number: 


k = — i—n, n 6 Z. (18) 

o 

I.e. for k 7 ^ 0 there is an instability increasing either to the right or left relative to the 
flow. Nonzero modes will not be excited if 6 is sufficiently large. We can compare 5 with 
the parameters of liquid in three ways, but since this instability is due to hydrodynamic and 
electro dynamic forces then we choose the option for 6 where v and v m are included equally. 
Therefore, considering the convection assumption we have a range for 5: 


sJATXpVV n 


«( 5 < 


y/Anpu n 


B 0 B 0 

which implies that v v m or magnetic Prandtl number Pr m -C 1 and also Ha 1. 


(19) 


4 











One-dimensional perturbations. 


It should be noted that for k 7 ^ 0 the obtained system can be solved exactly, but we restrict 
our investigation to a one-dimensional flow in the y-stable region. Then vector perturbations 
are two-dimensional (no z component). We obtain two independent systems for potentials 
( 1201 ) , (EH) and for the components of perturbations Q 22 ]) , (1231) : 


x'ya 


iiq - 


i^v — 


the pressure disturbance is expressed through disturbance 


R dq 
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( 20 ) 

Bq da 

d 2 q 

= 0 , 

( 21 ) 
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( 22 ) 
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(23) 

47 t p dz 

V dz 2 




Bb 

^ = ~47T 

Boundary conditions are as follows: 

M±<5) = V(±S) = ||(±<5) = ^(±5) = 0. (25) 

Note that when v = v m there is a symmetry in both systems: after a replacement b —y 
V\JAnp (or a —y q^JAnp) the form of both equations doesn’t change. Therefore, the resulting 
spectrum is degenerate, so the eigenfunctions can be found in the form b = av^/Airp, where 
a = ±1. We obtain the eigenvalues corresponding to a stable flow: 


(24) 


il = -77TK?! 71 ' 2 ™ 2 + H q2 ) < 0) n G Z. ( 26 ) 

(2 sy 

This means that for the close values of kinematic and magnetic viscosity small one¬ 
dimensional perturbations are damped and the flow is stable. However, as follows from ffl9]) . 
there already exist some undamped perturbation modes. 

Solution of the system (l 20 lh (| 2 B and (| 22 ll . (1231) gives us the same eigenvalues A and leads 
to the same equation for the eigenvalues of the increment 7 . We solve (l22]h (1231) by seeking 
a solution in the form b = boe Xz , v = i>oe Xz . It leads to the eigenvalues: 
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2 
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Then we have a solution in the form: 


b(z) = bi(C ie x " + C 2 e- A ”) + b 2 (C 3 e A « + C 4 e- Aa '), 


(27) 


(28) 
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Figure 2: The real roots of the equation (135|) . Solid line: i'yy/i'y + Ha 2 . Dashed line: 
tanh a /?7 + Ha 2 . 


v{z) = vi(C\e XlZ - C 2 e~^ z ) + v 2 (C 3 e A2Z - C 4 e~ X2Z ), (29) 

where the eigenvectors are chosen in the following way: b\ = b 2 = Bq\ 2 8 : V\ = 

5(27 — z/ m Ai 2 ), v 2 = 5(27 — r , m A 2 2 ) (Ai taken with the positive sign before the square root). 
With such a choice the coefficients C are dimensionless. 

Due to the boundary conditions we obtain a homogeneous system of four linear algebraic 
equations, which admits nontrivial solutions for A satisfying the equation: 

b\v 2 cosh(Ai<5) sinh(A2<5) = v\b 2 sinh(Ai<5) cosh(A 2 5). (30) 

Equation (|30l) can be factorized: 


b 4 v 2 = B 0 Vo tanh(Ai<5), (31) 

Vib 2 = B q Vq tanh(A 2 5), (32) 

cosh(Ai<5) cosh(A 2 <5) = 0, (33) 

where Vq is a quantity having the dimension of velocity. Choose for example Vq = 1 s/5, 


then 7 can be easily made dimensionless y5 2 /z/ —> 7. Eigenvalues A must satisfy either both 
equations (l3Tlh (1321) or equation (l33lh 

Because we restricted ourselves to the area where Pr m <C 1 we can make a Taylor expan¬ 
sion in terms of this parameter A. As a first approximation we set Pr m = 0. Then we obtain 
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the following eigenvalues: Ai<5 ~ y/i'y + Ha 2 , A 2 — 0. Equation (1331) gives stable roots: 

iq = —Ha 2 — + 7 m) 2 <0, n G Z. (34) 

The system (15TD . (1521) becomes a single equation, because the latter is satisfied identically: 

*7\/*7 + Ha 2 = tanh \j'i'y + Ha 2 . (35) 

One of the roots (stable one) is immediately visible: i'y = —Ha 2 < 0. Another root 
corresponds a pure imaginary increment indicating the flow instability: i'y > 0 (see Fig. [2]). 
Numerical solution of the equation (|35li gives us other roots among which there are those 
corresponding to the instability. 


Discussion and conclusions 

Fig. [3] shows a summary of the experimental results on the study of the stability of 
the Hartmann flow [5j. In these experiments the resistance coefficient A = —2 p'5/pV 2 has 
been measured. For the Hartmann flow it has the form: A h — 2 Ha/Re. At the figure the 
deviation from a bisector of the coordinate angle means that the flow for the given parameters 
is already turbulent. 



0 2.5 5.0 7.5 10.0 12.5 15.0 


Ah ■ IQ 3 


Figure 3: Experimental data on the resistance coefficient in comparison with the theory 
for Hartmann flow. Highlighted region denotes the range of applicability of the convective 
approximation. 

From (J4]) follows that Ha\/Pr m <C 1, so in the convective approximation A <7 Ha~ l Re m 
Since the experiments were carried out at Ha ~ 10 3 and Re m ~ 10” 1 (mercury) it is possible 
to say that A h ^ 10 -2 . And as can be seen from Fig. [3] in this held the how is not laminar. 
Thus the found instability really takes place. 

Fulfillment of the condition (1T91) provides only the absence of perturbations aimed parallel 
to the planes (i.e. the absence of x and y components of perturbation). As can be seen from 
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(T35l) with an increase of Ha the instability becomes suppressed. This can be explained by the 
fact that with increasing magnetic field (and as a consequence an increase in the Hartmann 
number) embeddedness effect begins to dominate, making it difficult to form instabilities 
in the initially laminar flow. But the equation (|35]) is obtained in the limit Pr m —> 0. 
To determine more precise conditions under which the flow could be sustained one needs a 
detailed analysis of ()30lh 

It should also be noted that in the limit Ha 0 the instability is still present, although 
the flow itself becomes a Poiseuille flow. Absence of such a limit transition was noticed 
in [2] and there this is explained by the fact that a non-zero magnetic field already generates 
instabilities which then can develop without the participation of the magnetic field. 
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